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ABSTRACT 

While the low dispersion IUE NEWSIPS data products represent a significant im- 
provement over original IUE SIPS data, they still contain serious systematic effects 
which compromise their utility for certain applications. We show that NEWSIPS low 
resolution data are internally consistent to only 10-15% at best, with the majority of 
the problem due to time dependent systematic effects. In addition, the NEWSIPS flux 
calibration is shown to be inconsistent by nearly 10%. 

We examine the origin of these problems and proceed to formulate and apply algo- 
rithms to correct them to ~ 3% level - a factor of 5 improvement in accuracy. Because 
of the temporal systematics, transforming the corrected data to the IUE flux calibration 
becomes ambiguous. Therefore, we elect to transform the corrected data onto the HST 
FOS system. This system is far more self-consistent, and transforming the IUE data to 
it places data from both telescopes on a single system. 

Finally, we argue that much of the remaining 3% systematic effects in the corrected 
data is traceable to problems with the NEWSIPS intensity transformation function 
(ITF). The accuracy could probably be doubled by rederiving the ITF. 

Subject headings: instrumentation: spectrographs — methods: data analysis — tech- 
niques: photometric — ultraviolet: general 
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1. Introduction 

This paper addresses problems with the absolute 
flux calibrations, thermal corrections, and time-dependent 
sensitivity corrections of the IUE (Boggess et al. 
1978) low-dispersion NEWSIPS ("final archive") data 
(Nichols & Linsky 1996, NL hereafter). We demon- 
strate that there are systematic errors of up to 15 % in 
NEWSIPS fluxes and describe correction procedures 
which reduce the systematics to a level compatible 
with the best possible signal-to-noise ratio achievable 
by a single IUE spectrum, i.e, ~ 3%. 

While carrying out an independent research pro- 
gram which involved fitting Kurucz (1991) model at- 
mospheres to NEWSIPS low-dispersion data of B- 
type stars (Fitzpatrick & Massa 1999), it became ap- 
parent that the absolute flux calibration was suspect. 
The basis of our suspicion was that the model fit resid- 
uals were large, strongly wavelength dependent, and 
independent of stellar spectral type. Since we were 
fitting energy distributions of main sequence stars 
throughout the range 10,000 <T eff < 30,000 K, it 
was difficult to imagine a single ionic signature that 
would produce such an effect. Consequently, we per- 
formed a detailed assessment of the NEWSIPS low- 
dispersion data, ultimately involving more than 4600 
spectra. This investigation revealed that not only is 
the absolute flux calibration of the NEWSIPS data 
inconsistent with its proposed standard, but that the 
data also contain serious thermal and time-dependent 
systematics. 

The nature of the NEWSIPS absolute flux calibra- 
tion problem is illustrated in Figure [l], which com- 
pares the available IUE NEWSIPS SWP and LWP 
data (see §|J) , HST FOS data, and models for the hot 
white dwarf G191B2B. The top spectrum is a ratio of 
the IUE NEWSIPS data and the Finley & Koester 
(FK) model used to define the IUE calibration (NL). 
There is obvious disagreement between the model and 
the NEWSIPS data, averaging ~ 5% and as large 
as 10% in the LWP. Furthermore, the wavelength- 
dependent structure in the ratio is nearly identical to 
the T e //-independent residuals we observed in our B 
star fits — including the high frequency "noise" vis- 
ible in the figure — indicating that this results from 
systematic calibration errors and not from random 
noise. The middle plot is the ratio of the FOS fluxes of 
G191B2B divided by the model provided by Koester 
to calibrate the FOS (Bohlin 1996, B96 hereafter). In 
this case, the agreement is excellent, but not exact, 



since G191B2B was only one of 8 white dwarfs used 
to calibrate the FOS (B96). The bottom plot is a ra- 
tio of the NEWSIPS and FOS data. It is very similar 
to the top ratio, but is also displaced downward by a 
wavelength independent, gray, displacement of ~ 2- 
3%. The displacement results from the IUE and FOS 
projects adopting slightly different models and differ- 
ent scalings for their UV calibrations; FOS is based 
on optical photometry calibrations (B96), and IUE 
on UV calibrations (NL). 

An example of the systematic time-dependent ef- 
fects present in NEWSIPS data is illustrated in Figure 
U for the IUE LWR camera. This plot shows the mean 
LWR flux in the wavelength band 2400 < A < 2800 
A for 3 standard stars as a function of time during 
the years when the LWR was the default long wave- 
length camera. The data for each star were normal- 
ized by the overall mean for that star during the en- 
tire time interval. Notice that each star shows the 
same time-dependence, indicating that the effect is 
instrumental, and not intrinsic to a particular ob- 
ject. The magnitude of the effect over this period 
is roughly 10%. The LWR camera shows the largest 
level of time-dependence. A preliminary discussion of 
the time dependent behavior of both the NEWSIPS 
low and high dispersion data can be found in Massa 
et al. (1998). 

High frequency ( "noise-like" ) temporal systematics 
are also present in the data. Figure || illustrates these 
effects. We first normalize the large aperture spectra 
of 3 IUE standard stars (HD 60753, BD+28°4211 and 
BD+75°375 - see §3) by their mean spectra. The nor- 
malized spectra were then arranged chronologically. 
The first 40 spectra of this set were divided into 2 
20-spectrum samples - each composed of every other 
spectrum. This insures that the 20-spectrum sam- 
ples have effectively identical temporal distributions, 
so time dependent systematics affect their means in 
the same way. The top curve in each panel is a ra- 
tio of the 2 20-spectrum means. The middle curve in 
each panel is a ratio of 2 20-spectrum means obtained 
from the last 40 spectra of the same 3 star sample and 
similarly prepared. Finally, the bottom plot in each 
panel is a ratio of the 2 40-spectrum means obtained 
from all of the spectra used in the first and second 
curves. While time dependent effects will cancel in 
the first 2 curves, they are maximized in the last one, 
since it compares data obtained at the beginning and 
end of the mission. If the high frequency structure 
was due to random noise, its amplitude in the bottom 
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curve would be 1/ v2 times smaller than in the top 
two curves, since twice as many spectra are used to 
produce the ratio. However, the amplitude of the high 
frequency structure is clearly larger in the bottom ra- 
tio! This is because much of the structure is not due 
to random noise, but to time dependent effects in the 
data. Figure || demonstrates that the peak-to-peak 
amplitude of the high frequency systematics often ex- 
ceeds 10%. 

In addition to temporal effects, residual tempera- 
ture effects are also present in the NEWSIPS data. 
These are characterized by the camera head amplifier 
temperature, THDA, measured at the beginning of 
the exposure (these values are supplied as part of the 
NEWSIPS headers, see, Garhart et al., 1997). The 
THDA effects are localized in wavelength. To demon- 
strate the effect, we used the same data set described 
in the previous paragraph, and binned them over a 
small wavelength band. These data are then plotted 
in chronological order in the top portions of Figures 
and b. The same data are then rearranged by 
THDA value and plotted in the lower portions of the 
figures. The presence of a systematic THDA effect on 
the order of 5% is obvious. We only show the LWP 
and SWP data, since the large temporal effects in the 
LWR data tend to obscure the smaller THDA system- 
atics. 

The previous examples demonstrate that the IUE 
absolute flux calibration is inconsistent with its refer- 
ence model by as much as 10%, that the LWR data 
contain time-dependent errors of similar magnitudes, 
that all of the data contain high frequency tempo- 
ral effects whose amplitudes exceed 10%, and that 
thermally induced systematics on the order of 5% are 
also present. We emphasize that systematic errors of 
this order can be important in many applications and 
must be corrected for the following reasons: 

1. To utilize the unprecedented temporal baseline 
of IUE, time dependent instrumental drifts in 
the data must be corrected to the level of the 
maximum achievable S/N. 

2. To take full advantage of the signal-to- noise 
(S/N) capabilities of IUE, an absolute calibra- 
tion whose uncertainty is equal to or less than 
the maximum achievable S/N is required. 

3. Time dependent dependent high frequency struc- 
ture can mimic the strengthening or weaken- 
ing of spectral features. These effects can pro- 



duce misleading interpretations and act to nul- 
lify potential noise reductions gained by averag- 
ing spectra. They must be reduced to the level 
of typical point-to-point noise. 

So what is the maximum achievable S/N for IUE 
data? In spite of the limited dynamic range of the 
IUE detectors, we demonstrate in §3 that a S/N ~ 
30 : 1 should, in principle, be possible for a single 
spectral resolution element in an optimally exposed 
spectrum. Furthermore, many bright objects were 
either observed more than once or else observed in the 
trail mode, which could increase the S/N by a factor 
2 (see §3). As a result, a large fraction of the objects 
in the IUE archive have a potential S/N ~ 30 : 1 - 
roughly 4-5 times more accurate than the current 
NEWSIPS calibrations. 

In the following sections we examine more closely 
the systematic effects present in IUE NEWSIPS data 
and present a scheme for the removal of these effects 
and a correction to the absolute calibration. Bless 
and Percival (1998) performed a critical review of the 
available UV calibrations and deduced that the FOS 
absolute calibration is superior. Consequently, we 
elected to derive a transformation between the IUE 
NEWSIPS and FOS systems, rather than recalibrat- 
ing IUE using the G191B2B model. This also insures 
that both datasets are on a common scale. 

In §|^, we introduce the stars to be used in our 
analysis and indicate which ones will be used to de- 
termine temporal and thermal trends, to derive an 
absolute flux calibration, and to verify the results. In 
§||, we describe some basic characteristics of IUE and 
explain the different observing modes and why each 
must be calibrated separately. We also discuss how 
the available data were culled into our final sample. 
In §Q, we describe the mathematical formulation of 
how we correct for systematic effects. In ^5|we apply 
the analysis to the IUE data and present our results. 
In we verify our results by applying them to se- 
quences of spectra not included in the derivation of 
the corrections. In we analyze both the random 
and systematic errors present in the corrected data. 
In §||, we summarize our conclusions and discuss the 
availability of IDL programs which apply the correc- 
tions to NEWSIPS low dispersion data. 

2. The program stars 

Basic data for the program stars used in this study 
are listed in Table 1. The stars were selected from the 
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IUE standard stars compiled by Perez et al. (1990) 
and from the FOS calibration stars given in B96. 
We also indicate in the table the role of each star 
in our analysis. The "temporal/thermal standards" 
are used derive the time- and temperature- depen- 
dence of the instrumental response, and the "tempo- 
ral/thermal control" stars are used to verify the re- 
sults. The "flux standards" are the stars used to de- 
rive the transformation between the IUE NEWSIPS 
and FOS flux scales and the "flux control" stars are 
used to verify these results. 

3. The data 

In this section, we describe the selection of the 
spectra used in our analysis. We begin with a discus- 
sion of a few general characteristics of IUE data and 
observations (§3.1) and then describe the 4 IUE ob- 
serving modes (§3.2). Lastly, we list the criteria used 
to select or reject individual spectra and individual 
data points (§3.3). For a more detailed description of 
the general properties of IUE and its data, see New- 
mark ct al. (1992) and Garhart ct al. (1997). 

3.1. Overview 

IUE had two UV spectrographs, covering the short 
and long wavelengths regions. Each spectrograph 
could send its output to either of 2 (primary and re- 
dundant) cameras. Consequently, IUE spectra are 
referred to as long wavelength primary (LWP), long 
wavelength redundant (LWR) or short wavelength 
prime (SWP). The short wavelength redundant cam- 
era never operated properly and there are no NEWSIPS 
data for it. The wavelength coverages of the spectro- 
graph/camera combinations were 1150 < A < 1975 A 
for the SWP, 1910 < A < 3300 A for the LWP and 
1860 < A < 3300 A for the LWR. 

Each spectrograph could be accessed through one 
of 2 object apertures. One was a 10 x 20 arc sec oval, 
the large aperture, and the other was a 3 arc sec circle, 
the small aperture. The size of the telescope image in 
the aperture plane was ~ 4 arc sec at best, and thus 
overfilled the small aperture by a large margin (see 
Garhart ct al. for a more detailed discussion). 

The IUE detectors had 768 x 768 pixels and each 
pixel had a very limited dynamic range, with only 
256 distinct output levels (8 bits). Since IUE was in 
high earth orbit, part of the output signal was due to 
particle background, as well as to a pedestal of read 
noise. Consequently, the maximum S/N possible for 



a single pixel was ~ 10 : 1 (see Ayers 1993 for a more 
thorough discussion). Since a low dispersion spec- 
trum typically had a full width perpendicular to the 
dispersion of about 3 pixels, the maximum S/N of a 
spectral pixel was ~ \/3 x 10 or ~ 17: 1. Furthermore, 
<~ 3 spectral pixels make up a spectral resolution el- 
ement, so the maximum S/N for a single spectral el- 
ement was ~ V3 x V3 x 10, or ~ 30 : 1. However, 
such high S/N rarely occurred over a large spectral 
range since the limited dynamic range of the detec- 
tor made it necessary to under expose some pixels in 
order to avoid saturating others. Nevertheless, when 
multiple exposures of the same object are averaged, a 
S/N ~ 30 : 1 should be attainable over broad spectral 
regions. Therefore, our goal is to reduce the system- 
atic effects in the IUE data to a similar level, i.e., 
-3%. 

Once an exposure was obtained, the IUE cameras 
had to be read and then prepared for the next expo- 
sure. This "read-prep" operation typically required 
about 20 minutes. This was often far longer than the 
actual exposure times and several of the observing 
strategies described below were intended to minimize 
the impact of the read-prep overhead time. 

Finally, the IUE satellite locked onto its target us- 
ing its Fine Error Sensor (FES). In this mode, the 
satellite could point to a fixed position on the sky 
with an accuracy of 0.25 arc sec. 

3.2. The observing modes 

IUE spectra were obtained in one of 4 observing 
modes, each of which had its own unique advantages. 
These 4 observing modes were: 

Large Aperture: This mode simply involved ob- 
taining spectra by centering the object in the large 
aperture and was the primary observing mode. Since 
the large aperture was ~ 3 times larger than the im- 
age of a point source produced by the IUE optics, the 
0.25" tracking capability of IUE was not critical for 
large aperture exposures, and the photometric quality 
of these data is the best. 

Small Aperture: In this mode, the star was cen- 
tered in the small aperture. Since the image was con- 
siderably larger than the aperture, the exact amount 
of flux entering the detector was highly dependent 
upon the pointing. Consequently, the overall pho- 
tometric quality of these data is very poor, varying 
by as much as a factor of two. Nevertheless, small 
aperture could be valuable for several reasons. First, 
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low dispersion large and small aperture spectra were 
well-separated on the detector and so both could be 
recorded without an intervening read-prep. Thus the 
overhead requirements and poor dynamic range of the 
detectors could be circumvented to some degree by 
exposing the large and small aperture spectra to dif- 
ferent signal levels and later combining the extracted 
data. (Small aperture data would be scaled to the 
photometric level defined by large aperture using data 
in mutually well-exposed wavelength regions.) Sec- 
ond, small aperture spectra had slightly better spec- 
tral resolution, due to the smaller image on the de- 
tector. Third, sometimes it was necessary to use the 
small aperture in order to isolate objects in crowded 
fields. 

Trailed: In this mode, the star was allowed to drift 
across the large aperture, in the cross-dispersion di- 
rection, during an exposure. This produced a widened 
spectrum which had 2 advantages. First, since nearly 
4 times as many pixels were exposed, trailed spectra 
have nearly twice the S/N of large aperture spec- 
tra with the same read-prep time. Second, since 
more pixels contribute to a single resolution element, 
trailed spectra are less sensitive to localized detec- 
tor irregularities (fixed pattern noise). On the other 
hand, the photometric accuracy of trailed spectra is 
somewhat inferior to that of large aperture spectra 
because the exposure time depends on the exact tra- 
jectory of the object through the aperture and the 
exact drift rate. Further, during the maneuver, the 
spacecraft had to rely on its gyros for stability and 
could not take advantage of the FES feedback. As 
a result of these effects, trailed spectra typically have 
the best relative photometric accuracy, but their over- 
all flux level is slightly less accurate than large aper- 
ture spectra. 

Multiple exposures: In this mode, the star was 
placed at 2 or 3 distinct locations perpendicular to 
the dispersion in the large aperture. The net result is 
a widened spectrum, although it is not uniformly ex- 
posed (hence, this mode was often referred to as the 
pseudo-trail mode). Since FES tracking was in effect 
throughout the exposures, the photometric accuracy 
of these spectra are probably better than for trailed 
spectra, but the standard stars were not observed of- 
ten enough in this mode to verify this supposition. 

For the purposes of this paper, 2 important points 
emerge from the preceding discussion. 

1. Each observing mode exposed different portions 



of the detector. Therefore, the temporal and 
thermal behavior and absolute calibration of 
each mode must be considered separately. Un- 
fortunately, we lack the data to do this for the 
pseudo-trail mode and it will have to be cali- 
brated from the other modes. 

2. Errors in IUE spectra typically contain two dis- 
tinct components: point-to-point (or relative) 
errors which are important for measuring spec- 
tral features, and scaling errors, which affect the 
overall level of the spectra and are important 
in fitting models or in concatenating IUE spec- 
tra with each other or with optical photometry. 
The relative magnitude of these 2 types of er- 
ror depends upon the observing mode. Further- 
more, the scaling errors (which originate from 
pointing and focus inaccuracies) can sometimes 
be quite large. 

There are additional effects which can have a 
strong influence on IUE spectra. For instance, the 
particle background rate could sometimes become 
rather large, reducing the S/N of spectra of even 
bright objects to only a few. There were also non- 
random effects which included microphonics noise (re- 
ferred to as "PINGS") and minor frame telemetry 
drop outs, which typically have a strong affect on a 
localized portion of the spectrum. See Garhart et al. 
(1997) for a more complete discussion of these and 
other effects. 

3.3. Data selection 

Table 2 contains information on the IUE NEWSIPS 
data used for this analysis. We submitted a request 
to the NSSDC for all of the available NEWSIPS low 
resolution spectra for the stars listed in Table 1 at 
the end of 1998 March. All but a few spectra were 
available at that time. The delivered data were then 
screened as follows: 

1. The ranges of exposure times were restricted de- 
pending on the object, camera, observing mode, 
and application. In Table 2, t m i n and t max give 
the lower and upper limits, respectively, on ac- 
ceptable exposure times. 

2. The LWR sample was restricted to the time pe- 
riod when the LWR was the default long wave- 
length camera (1978 - 1984). 
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3. Outliers were rejected. These were defined as 
follows: The mean flux for each spectrum was 
determined over a pre-specified wavelength in- 
terval (1400 - 1700 A for the SWP and 2400 - 
2800 A for the LWR and LWP). The sample 
mean and RMS scatter (a) were determined for 
these mean fluxes. If the absolute value of the 
difference between the mean flux for a given 
spectrum and the sample mean differed by more 
than 3cr, then the spectrum was rejected. This 
criterion was applied iteratively until no addi- 
tional spectra were rejected. 

4. A few LWRs were rejected "by hand" because, 
although they passed the outliers criteria, their 
shapes were distinctly peculiar. All of these 
were observations of HD 60753 and most of 
them were affected by "pings" (data drop outs) 
or extreme background levels. 

The number of spectra which survived the screen- 
ing process outlined above is listed in the last column 
of Table 2 for each star and observing mode. Within 
each acceptable spectrum, data points with i^-Flag 
values (see Garhart et al. 1997) equal to (no known 
problem), —128 (on the positively extrapolated ITF) 
or —256 (on the negatively extrapolated ITF) were 
given weights of unity and all other points were as- 
signed zero weight. 

Finally, we note that there are two distinct ITFs 
for LWR spectra (see Garhart et al. 1997) and these 
result in slightly different wavelength scales. Fur- 
ther, ITF A (identified as LWR83R94 in the MXLO 
headers) results in spectra with 563 data points while 
ITF B (LWR83R96 in the headers) spectra have 562 
points. We used the wavelength scale from ITF A 
for all of the LWR spectra. This ignores a difference 
between the 2 wavelength scales which increases lin- 
early from 0.2 A at 1950 A to 1.66 A at 3150 A (the 
longest wavelength we calibrate). However, since even 
the largest deviation is smaller than the sampling in- 
terval (2.67 A) and much smaller than a resolution 
element (~ 7 A at the longer wavelengths), we felt 
that interpolating the data from one grid to the other 
was unwarranted. Therefore, we simply accept the 
minor systematic error which arises from adopting a 
common wavelength scale for all of the spectra. 



4. Mathematical description of the analysis 

To bring the NEWSIPS data onto a common scale 
with the FOS data, we must first remove its time and 
THDA dependencies. We do this by fitting the depen- 
dencies and then applying the results to correct the 
data for the systematics. We then derive the trans- 
formation between the corrected data and the FOS 
system. In this section we provide a mathematical 
outline of the problem. We describe its application to 
the data in §||. 

4.1. The corrections 

We wish to analyze {i — 1, . . . , M} standard stars 
observed at different wavelengths and times to deter- 
mine the time degradation and THDA dependence of 
the instrumental response. In the analysis, we adopt 
the following model of the temporal, t, and THDA, T, 
dependence for the i th standard: 

/(A, t, T)i = /(A, t , T a ) t g(X, t - t )h(X, T - T ) (1) 

where g{X,t — to) and h{X,T,To) are assumed to be 
universal multiplicative functions which describe the 
time and THDA dependencies of the instrumental re- 
sponse at A and are equal to 1 at (i, T) = (to, To). 
The fact that we have written the t and T dependen- 
cies as separate functions implicitly assumes that the 
form of the THDA dependence does not change with 
time and that the temporal dependence is the same 
for all values of THDA. 

Taking logarithms linearizes the problem, viz., 

log /(A, t, T)i = log /(A, t , T ) 4 + log <?(A, t-t ) + 

\ogh(T~T ) (2) 

A simple form for the functions log g(X,t — to) and 
logft,(A,T — To) which satisfy our assumptions is 

K 

\ogg(X,t-to)=Y,a( X )k(t-to) k (3) 

k=l 
L 

\ogh(X,T-T ) = Y,bWk(T-To) 1 (4) 
l=i 

i.e., K and L th order polynomials. 

It is possible to fit the data set for each stan- 
dard star individually using the previous equation, 
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thereby obtaining the two sets of coefficients, {a(X)k} 
and {b(X)i}, and the flux at the fiducial values, 
log/(A, to)i, for each standard star separately. How- 
ever, we would like to fit the data of the M stan- 
dards simultaneously, thereby determining a univer- 
sal estimate of the coefficients, utilizing all of the 
available data. To accomplish this, we first concate- 
nate the data into a single data set. If there arc 
{rn = 1, . . . , M} standard stars, each with N m obser- 
vations at times {t n \n = 1, • • • , N m }, at each wave- 
length A, then the concatenated series {y(X, t)} is de- 
fined as, 



{y(X, t, T)} = {log /(A, ti)i, . . . , log /(A, t Nl )i, 
log/(A,ti) 2) • • ■ ,log/(A,tjv 2 )2, 
• • .,logf(X,ti) K , 

logf(X,t NM ) M ] (5) 

The temporal and THDA dependence of the combined 
data set at each wavelength is then fit with a standard 
linear regression model of the form 



y(X,t,T) = Epilog /(A, t ) m X 



Om 



Ek=i*WkXk + ZLibWiYi (6) 

where the Xo m are "box-car" functions which are ei- 
ther or 1, depending on whether the data refer to 
the m th star, the X/. are polynomials of the form 
[t — to) k , and the are polynomials of the form 
(T - T ) k , where the t and T are the time and THDA 
corresponding to the particular term. 

As pointed out in §||, the major source of error 
in the spectra is often an overall scaling factor due 
to inexact centering of the object in the aperture or 
slight trailing errors. To suppress this effect, we nor- 
malized the spectra by their mean flux over a wave- 
length band Ai < A < A2. These normalized spectra 
are denoted as r(A, t, T). The r(A, t, T) are fit at each 
wavelength by equation ^, and then the normalization 
constants arc fit independently the same way. As a 
result, we determine 3 sets of coefficients, the {a{X)k} 
and {6(A);} in equation || (except now they apply to 
the r(A, t, T)) and a set {a k',b i}, which fit the level 
of the flux in the standard band, relative to its value 
at t — to, T — T . Consequently, to correct the flux 
of an object observed at time t with THDA = T, one 
must divide the observed flux by the function: 



K 

g(X,t-t )h(X,-T ) = f[ 10 [Q(A)fc+aofcl(t ~* o) x 

k=l 

Y\ lOll>Wi+boi](T-To) 1 Pj\ 
1=1 

The result is how the spectrum would have appeared 
if it had been obtained at t = to with THDA = T . 

Finally, since there are not enough pseudo-trail 
spectra to perform an independent calibrat ion , this 
case is treated differently, and discussed in §5.2. 



4.2. Flux scale transformation 

Once the temporal and thermal corrections are de- 
termined, the transformation to the FOS flux scale is 
relatively straightforward. It is simply a mean of the 
ratios of the FOS spectra of the "flux standards" (re- 
binned to IUE resolution) to their mean IUE spectra. 

5. Application of the analysis 

In this section we provide the details of the general 
analysis outlined in the previous section, as applied 
to the data for the program stars listed in Table 1. 

5.1. Temporal and THDA corrections 

In performing the fits, we used sixth degree poly- 
nomials in both t and T for the g(X,t) and h{X,T). 
The time, t, was expressed in Julian years, and to = 
2445000/365.25 (1 Feb. 1 1982). The fiducial THDA 
value, To, was set to 9 for the SWP and LWP data 
and 13 for the LWR. Due to a paucity of data at ex- 
treme THDA values, THDAs of SWP and LWP data 
less than T m i n = 6 were set equal to 6 and values 
greater than T max = 13 were set equal to 13. The 
same process was used for the LWR data except with 
T m in = 11 and T max = 16. 

For normalizing the spectra to obtain r(A, t, T), we 
used 150 data points in the range 2399.69 < A < 
2796.43 A for the LWP and LWR and 179 points in 
the range 1400.35 < A < 1698.74 A for the SWP. 

The stars HD 60753, BD+28°4211 and BD+75°375 
were used as primary standards for the corrections 
because they have the largest number of spectra, and 
these span the entire lifetime of IUE. BD+33°2642 
and HD 93521 have the next largest number of spec- 
tra (roughly half of any one of the primary standards) , 
and excellent temporal coverage. These objects will 
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be used to verify the results derived from the 3 stan- 
dards (see §0). 



Figures 5a -c give examples of the fits to the time 
trends in the relative scale factors and at the specific 
wavelengths for the 3 IUE cameras. The solid curves 
are the fits for THDA = To, and the dashed curves 
are for the cases T - T Q = 0.75{T max — T ) and T - 
T = 0.75(T min - T ). The fits for the extreme THDA 
values usually parallel the T = To curve and typically 
represent a much smaller effect (as expected from our 
discussion in §[!]). Data from the 3 standards used to 
determine the fits are depicted by different symbols. 
Each of the individual wavelength plots shown are 
actually the means of 3 adjoining wavelength points, 
to reduce the overall noise. It is clear that the data 
for the 3 stars arc completely interspersed and that 
the solution is consistent with all three. Several other 
aspects of the plots are of interest. 

1. There are no major time systematics in the 
LWP data (Fig. §). 

2. Data longward of 3100 A become very unre- 
liable in both of the long wavelength cameras 
(Figs.^-b). 

3. The signal-to-noise of data shortward of 2000 
A is very poor in the LWP data, but relatively 
good in the LWR (Figs. |a|-b). 

4. The large time dependent systematic in the scal- 
ing of the LWR data is clearly demonstrated 
(Fig. H). 

5. All but the shortest wavelengths of the SWP 
data have comparable signal-to-noise (Fig. |k]). 

6. Strong systematic effects are present in the 
short wavelength SWP data (Fig. |5c| ). 

7. The SWP data are generally of higher quality 
(recall that these are all comparable exposures). 
Both the scale factors and the individual wave- 
length fits have smaller dispersions in the SWP 
data. 

5.2. Flux transformations 

Once the spectra are corrected to their fiducial 
time and THDA values, the transformation to the 
FOS system is straightforward. The stars used to de- 
termine the flux transformation were BD+28°4211, 
BD+75°375, BD+33°2642 and G191B2B. Both FOS 



and high quality IUE data are available for each of 
these. HD 60753 was not observed with FOS, but will 
provide a powerful verification of the flux calibration. 

The FOS data were first smoothed to the IUE res- 
olution using the data provided by Garhart et al. 
(1997). These spectra were then sampled onto the 
IUE grid. There is, however, one complication. In 
order to make the sharp He I features located through- 
out the long wavelength IUE spectra of BD+75°375 
cancel with their counterparts in the FOS spectra, it 
was necessary to adjust the wavelength scale of the 
long wavelength cameras. Since experience has given 
us considerably more confidence in the FOS calibra- 
tions, we adopted the FOS wavelength scale and de- 
rived a set of adjustments for the IUE scale. The 
measured differences are listed in Table ||. In prac- 
tice, we applied a spline interpolation between these 
points. 

Figures ^-c show the ratios of the completely 
corrected large aperture data (with the wavelength 
corrections applied to the long wavelength data) di- 
vided by the corresponding FOS data. Curves from 
the 4 primary standards, BD+28°4211, BD+75°375, 
BD+33°2642 and G191B2B are depicted by differ- 
ent line styles. The mean curve, used for the cali- 
bration, was formed by first adjusting all of the ra- 
tios to the sample mean value across the same wave- 
length bands described above and then determining a 
weighted mean ratio, where the weighting factors were 
just the number of observations that entered each ra- 
tio. The standard deviation of the weighted mean 
ratio was also calculated and is shown at the bottom 
of each plot. The excellent agreement of the differ- 
ent curves emphasizes the reality of the structure, in- 
cluding the large point-to-point structure. The RMS 
dispersion shows that the overall internal agreement 
of the calibration curves is ~ 1% - well within our 
goal. The trailed and small aperture ratios have sim- 
ilar scatters. It is interesting that the feature referred 
to as the 1515 A feature by Garhart et al. (1997) is 
clearly present. We shall return to this point later. 

As expected, the corrections are very similar to the 
curves shown in Figure |l|. In addition to a general 
gray offset of ~ 5%, there is also structure present at 
the ~ 10% level in all of the camera. 

5.3. Pseudo-trail spectra 

As mentioned in §[| the pseudo-trailed (p-trailed 
hereafter) spectra present a special problem since 
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there are not enough of them to perform a thorough 
analysis of their time and THDA systematics. Table 
| shows that there are only 18 LWP, 11 LWR and 35 
SWP p-trail spectra for the standards. These include 
spectra with both 2 and 3 exposures in the large aper- 
ture. These subsets do not expose the same pixels in 
exactly the same way and there is no a priori reason 
to assume their corrections will be similar. However, 
we are forced to assume that they are, since we lack 
the data to determine otherwise. This situation is un- 
fortunate, since in spite of the paucity of p-trail data 
for the standards, it was a popular observing mode, 
and there are many p-trail spectra in the archive. 

Due to the lack of data, we had to adopted the fol- 
lowing approach for calibration of the p-trail spectra. 
We assume that we know the intrinsic flux distribu- 
tions for the p-trail spectra from either FOS spectra 
(if available), or mean values of fully corrected IUE 
large aperture data, transformed to the FOS system. 
Each p-trail spectrum was then divided by its corre- 
sponding FOS or mean IUE large aperture spectrum 
to produce a set of normalized spectra whose mean 
value should be unity. We then corrected the normal- 
ized p-trail spectra with both the large aperture and 
trailed temporal and THDA corrections and IUE— to- 
FOS calibration and compared the results. 

We found that the large aperture corrections and 
calibration performed best in all cases; removing all 
obvious trends from the data, and reducing the over- 
all scatter. They also produced a mean which was 
uniformly close to unity. This result was somewhat 
surprising, since the p-trails expose a wide swath of 
pixels, and one might expect their properties to be 
more similar to trailed spectra. However, it is possi- 
ble that adjacency effects due to "beam pulling" dom- 
inate. If so, the multi-peaked cross-dispersion struc- 
ture of the p-trail spectra may make their properties 
most similar to a single large aperture spectrum. 

In any event, the RMS scatter of the normalized 
p-trail spectra (corrected by the large aperture rela- 
tionships) from unity is ~ 1% over most of the usable 
range of the SWP and LWP, but there are regions 
where systematic deviations of ~ 3% may be present. 
For the LWR, the the scatter is uniformly ~ 3%, but 
this is largely due to the overall poor photometric 
quality of the available LWR p-trails. 



5.4. Special wavelength regions 

There are wavelength regions for each camera where 
either the intrinsic data or our corrections algorithms 
are not well-defined. The wavelength extremes over 
which the corrections can be applied were determining 
by examining plots such as those shown in Figures 5a 
- [kj. These plots show that data for the longest wave- 
lengths of the long wavelength cameras are poorly 
defined, and applying correction factors to these data 
has little meaning. Table || lists the wavelength range 
over which our correction factors are considered re- 
liable for each camera. The factors are set to unity 
outside these regions. 

There are also some specific wavelength regions 
which are problematic. However, most of these should 
not be a real concern, since they are flagged by the 
j^-Flags as being poor quality data. However, if the 
correction algorithm is blindly applied, it will produce 
a number for such data. Therefore, we caution users 
of our correction scheme to always use the ^-Flags to 
eliminate problematic data points. 

One region which is not flagged by the i/-Flags, but 
is unreliable is the region near Ly a in the SWP data. 
The NEWSIPS spectral extraction uses a low order 
polynomial to represent the background on either side 
of the spectrum. This approach cannot handle Geo- 
coronal Ly a emission which fills the aperture in long 
exposures. As a result, long exposures will be contam- 
inated by Ly a emission over a spectral region equal 
to the projected size of the aperture used to obtain 
the data. These ranges are 1207 - 1222 A for large 
aperture are trailed data and 1210 - 1221 A for small 
aperture data. NEWSIPS low dispersion data cannot 
be trusted in these regions, and the corrections have 
been set to unity over them. 

Finally, we note that Garhart et al. (1997) demon- 
strate that the NEWSIPS SWP spectra show an 
anomaly near 1515 A. However, this problem appears 
to be significantly reduced in the corrected data, and 
we no longer consider it to be a major problem. 

6. Verification 

We now must verify the temporal and thermal cor- 
rections and the flux transformations derived in the 
previous section. In doing so, it is mandatory that 
we use only spectra that were not employed to derive 
the relationships. Since the best data sets were used 
to derive the relations, we cannot expect to test the 
full accuracy of the results. 
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We begin with verification of the temporal and 
THDA corrections. We use BD+33°2642 for verifica- 
tion since it and HD 93521 have the most data of the 
stars not used in deriving the temporal and THDA 
corrections. Figure |?] compares ratios of means of 
BD+33°2642 spectra obtained early in the mission 
to means of spectra taken late in the mission. Each 
set of ratios is labeled by the camera used to obtain 
the spectra. The dotted curves are ratios of uncor- 
rected NEWSIPS data and the solid curves are ratios 
of NEWSIPS data corrected for temporal and THDA 
systematics. The SWP and LWP ratios are 50 spec- 
tra means and the LWR ratios are 20 spectra means. 
The mean time of each spectral mean is provided on 
the ordinate labels. 

While (as expected from Figure || which uses far 
more data) the temporal corrections make little dif- 
ference in the LWP data, they have 2 effects on the 
LWR and SWP spectra. In each case, the ratios of the 
corrected data are closer to unity (much more so for 
the LWR data) and much of the point-to-point varia- 
tion is reduced in the corrected data, demonstrating 
that it was not true noise, but rather systematic ef- 
fects. 

We now turn to verification of the flux transforma- 
tion. Figures |Ia| - |8c| compare the fully corrected and 
uncorrected NEWSIPS data for all of our program 
stars. We begin by examining specific improvements 
in the spectra of stars used to derive the transforma- 
tions and then turn to those stars used to verify the 
results. 

The improvements for the stars used to derive the 
transformation (BD+75°375, BD+28°4211, BD+33 
2642, and G191B2B) are truly spectacular. In partic- 
ular, 

1. Reduction of point-to-point "noise" is most no- 
ticeable in the long wavelength cameras. In 
particular, the He I lines in BD+28°4211 and 
BD+75°375 are much more distinct. In fact, 
these lines are barely visible in the uncorrected 
NEWSIPS LWR spectra of BD+75°375, but are 
obvious in the corrected spectra. 

2. The 1515 A artifact (Garhart et al. 1997) 
is clearly present in the NEWSIPS spectra of 
BD+28°4211, BD+33°2642 and G191B2B, but 
it is reduced or completely removed in the cor- 
rected spectra. 



3. Structure in the region 2200 < A < 2500 A in 
the LWP spectra of G191B2B is removed. 

We must seek verification of these results in the 
stars which did not enter into the derivation of the 
relationships. We are at a bit of a disadvantage here, 
since the white dwarfs not included in the deriva- 
tions are not well observed and the OB stars have 
rather "busy" spectra. Nevertheless, the following are 
clearly seen: 

1. The point-to-point "noise" is clearly reduced in 
the corrected SWP spectra of HD 60753 and 
dramatically reduced in the LWP and LWR 
spectra of HD 93521 and HD 60753. 

2. The 1515 A feature clearly reduced in the SWP 
spectra of GD 153, GD 71 and HD 60753. It 
is also reduced in HD 93521, but it is difficult 
to see since it lies in a strongly blanketed re- 
gion of its spectrum. Given the vastly different 
flux levels and temporal distributions of these 
observations, the possibility that this artifact is 
completely removed by the corrections is quite 
good. 

3. The structure between 2200 and 2500 A in the 
LWP spectra of GD 153 and GD 71 is reduced 
in the corrected data. Its removal is not so ap- 
parent in LWP spectra of HZ 43 and HZ 44 be- 
cause of their higher noise level (see their FOS 
spectra). 

Finally, Figures |9"af|9c| compare the fully corrected 
NEWSIPS data with the FOS data. The FOS data 
have been degraded to match the IUE spectral resolu- 
tion. For the 3 calibration stars with the most data, it 
is almost impossible to distinguish between the FOS 
and IUE spectra. It is also clear that the corrected 
IUE spectra the 4 stars not used in the calibrations 
agree their FOS counterparts quite well. The only ex- 
ceptions are near Ly a in GD 153 and GD 71. That 
disagreement arises because these were relatively long 
exposures, so the region of Ly a is partly filled in by 
geo-coronal emission (see §5.4). 

We see, therefore, that the improvements provided 
by the new calibrations are also present in spectra 
which were not used to derive the relationships. This 
independent verification of our results provides confi- 
dence in their veracity. 
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7. Error analysis 

To quantify the significance of an observed feature 
or the accuracy of a flux level, two types of error must 
be evaluated: random and systematic errors. Broadly 
speaking, the random errors are due to uncontrollable 
effects which change in an unpredictable manner from 
one exposure to the next. They can be either point-to- 
point errors (e.g., photometric errors) or broad band 
errors (e.g., the scaling errors discussed in §|J). An 
important aspect of random errors is that they can 
be "averaged down", i.e., the average of N obser- 
vations repeated under similar conditions is 1/s/N 
times more accurate than any one of the observation. 
On the other hand, systematic errors depend upon 
some specific factor (e.g., exposure level), are typi- 
cally broad band, and cannot be averaged down, since 
the entire data set is subject to their influence. This 
does not mean, however, that they cannot be over 
come. The corrections we have derived for the tem- 
poral THDA systematics are examples of correctable 
systematic errors. 

7.1. random errors 

We begin by characterizing the random errors. 
This will first be done in a qualitative manner, us- 
ing the same approach adopted by NL. However, we 
use HD 93521 as our test object because it was not 
used in the derivation of either the temporal or flux 
corrections, making it an unbiased data set. Further, 
there are enough observations of the star to deter- 
mine whether the random S/N truly asymptotes. On 
the negative side, HD 93521 is known to have vari- 
able wind lines (e.g., Howarth and Reid 1993, Massa 
1995), and its spectrum is rather "busy", containing 
considerable structure. 

The S/N was calculated exactly as outlined by NL. 
First we summed n spectra (where n varies from 1 
to the the total number in the sample) drawn at ran- 
dom from the sample. Next, we calculated means and 
standard deviations over 4 point bins (roughly a reso- 
lution element). These were then converted into S/N 
ratios and summed over specified wavelength regions 
to obtain the final results. The wavelength regions se- 
lected were 1400 < A < 1500 A and 1650 < A < 1900 
A for the SWP and 2200 < A < 2900 A for the long 
wavelength cameras. These regions are comprised of 
the most responsive portions of the cameras and the 
SWP region avoids variable wind lines. 



the uncorrected and corrected NEWSIPS data. While 
there is only a modest improvement for the SWP 
data, the improvement for the long wavelength cam- 
era data is quit dramatic. The figure also shows that 
the LWR camera is the most "intrinsically" noisy of 
the 3, with a maximum attainable S/N ~ 40 : 1, fol- 
lowed by the SWP with a maximum S/N ~ 60 : 1 and 
the LWP being the best with a maximum S/N ~ 80: 
1. The figure also demonstrates that there is little to 
gain in summing more than ~ 10 IUE spectra, but 
up to that point, the gain is considerable. 

We next consider a more quantitative description 
of the errors. This is done by comparing the ob- 
served errors (the standard deviations derived from 
the repeated observations of the standard stars) to 
the NEWSIPS error model, whose results are given in 
the error vector in the NEWSIPS data files. Both of 
these are calculated as unweighted statistics, since the 
overall quality of the individual spectra are relatively 
uniform. Only data points without a known problem 
(f-Flag = 0) or on the extrapolated ITF (iv-Flag = 
-128 and -256) were included in each calculation. 



Figure 10 shows the results of the analysis for both 



Figure 1 1 shows ratios of the observed standard de- 
viations, a(Obs), to errors derived from the NEWSIPS 
error models, a(NEWSIPS), for the 3 standard stars; 
HD 60753, BD +28°4211 and BD+75°375. The 
a(NEWSIPS) are simply the square root of the quadratic 
mean of the NEWSIPS error arrays for all the good 
data points at each wavelength, while the a(Obs) are 
the standard deviations of all the good points at each 
wavelength. Each panel shows a different camera- 
observing mode combination. There are 2 curves in 
each panel. One is the mean a(Obs) for the 3 stan- 
dards using unsealed observations and the other is 
the mean a{Obs) derived from spectra which were 
rescaled to agree over the fixed wavelength bands de- 
scribed in §|^. For the large aperture data, these 2 
curves are nearly indistinguishable, since the scaling 
errors discussed in §|3] are nearly negligible. On the 
other hand, the two curves are well separated for the 
trailed data and very distinct for the small aperture 
data, with the unsealed observations always produc- 
ing larger errors. 

In every case, the NEWSIPS error model under- 
estimates the actual errors. The amount is typically 
of order unity for the large aperture and trailed data, 
but more than a factor of 2 for the small aperture 
scaled data. 

Table || lists the RMS scatter in the scaling fac- 
tor as a fraction of the flux across the bands given 
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in §|5| for each camera - observing mode combination. 
It also gives the mean a(Obs) /a(NEWSIPS) ratio for 
each camera - observing mode combination (with and 
without normalization) . These numbers will be useful 
guides when carrying out quantitative error analyses 
with NEWSIPS spectra, although it must be remem- 
bered that the ratios sometimes contain considerable 
shape, so their characterization as a single number 
can be an over simplification. 

7.2. systematic errors 

There are two additional parameters provided in 
with the NEWSIPS data files that can be used to 
search for systematic effects. These are exposure 
times and exposure levels. 

Since longer exposures typically have higher back- 
ground counts, it is possible that there could be sys- 
tematics in the data which are related to the exposure 
time used to obtain the spectrum. However, the re- 
sults of Figures |9a|-|9e] argue against such a systematic. 
The quality of the agreement between IUE and FOS 
data for the stars shown in the figure is excellent for 
stars with exposure times as short as 4 sec (for the 
SWP) and 6 sec (for the LWR and LWP) to as long 
as 21 minutes for the LWR and 30 minutes for the 
LWP and SWP. So we can be confident that the data 
are free of exposure time dependent systematics over 
this range of exposure times. 

Finally, we examined the data for systematic dif- 
ferences between spectra of the same star exposure 
to levels. For this purpose, we used the net spectra 
in the NEWSIPS data files, which are expressed in 
linearized flux units, FNs (see Garhart et al. 1997). 
Systematic differences between spectra exposed to dif- 
ferent mean FNs would indicate a problem with the 
Intensity Transfer Function (ITF) which transforms 
the observed counts (in data units, DNs), into the 
linearized FNs. 

To search for an ITF problem, we examined spec- 
tra of the same star obtained with different expo- 
sure times, making sure that saturated pixels were 
eliminated from the comparison. Although no major 
(> 3%) systematics were uncovered in the SWP and 
LWR data, the LWP spectra do contain sizable ITF 
systematics. Figure [l^ shows exposure level system- 
atics for LWP spectra of BD+28°4211, HD 60753, 
and BD+75°375. The plot shows mean NEWSIPS 
fluxes over the wavelength band 2350 < A < 2400A 
(the peak of the LWP camera response) for each star. 



These have been normalized by the mean of all expo- 
sures with 200 < FN < 400 (saturated pixels and pix- 
els using an extrapolated ITF were eliminated). Ob- 
servations of the same star with very different mean 
FNs result from different exposure times. The fact 
that fluxes derived from long exposure times (large 
FN values) are systematically larger (by about 5%) 
than fluxes derived from low FNs indicates a problem 
with the LWP ITF (since we have already ruled out 
systematics which depend solely on exposure time). 
This means that a comparison of an optimal expo- 
sure and a half optimal exposure of stars with similar 
energy distributions will contain systematic errors up 
~ 5%. 

Figure [l^ shows how the ITF problem can also af- 
fect the shape of an energy distribution. It displays 
ratios of long and short exposures of BD+28°4211 
(full curve), HD 60753 (dotted curve), and BD+75°375 
(dashed curve). The ratios consist of mean spec- 
tra derived from spectra whose mean FN over the 
band 2650 < A < 2700 A lie in the range 200 << 
FN >< 400, divided by means of spectra with 
600 << FN >< 800 over the same band. Data from 
saturated pixels and pixels using an extrapolated ITF 
were excluded from the means. It is clear from the 
figure how data near the peak of the camera response 
(between 2500 and 2850 A) are systematically differ- 
ent by up to 5%. 

The LWP ITF problem is not as severe as it first 
appears. It means that comparisons of well exposed 
and under exposed data may have systematics on the 
order of 5%. However, 5% is roughly the size of ran- 
dom errors for half optimal exposures (the reason the 
systematic effect shows up so well in Figure |l3| is that 
hundreds of spectra went into the ratios), so the effect 
does not dominate the errors when comparing a single 
half optimal exposure to an optimal one (although it 
is comparable to the random errors). Another prop- 
erty that tends to suppress the effect is that the shape 
of the response curve dominates the shape of the net 
spectrum. Thus, as long as the intrinsic shapes of 
the objects being compared are not too different and 
as long as the spectra have similar exposure levels, 
the impact of the LWP ITF problem should not be 
too bad. However, it underscores our assertion that 
the results of the present paper strictly apply to early 
type stars only. 

8. Summary and Conclusion 
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1. We have analyzed more than 4600 spectra to 
demonstrate that low dispersion NEWSIPS data 
contain systematic effects on the order of 10 — 
15% and to obtain corrections for these effects. 

2. Systematics were reduced to < 3% in most in- 
stances, but can be as large as 5% in a few spe- 
cific cases involving LWP data. Overall, we can 
hope for a S/N ~ 30 : 1 (~ 20 : 1 for some 
LWP applications) but not more. To exceed this 
value one would have to consider recalibrating 
the ITFs, and even rederiving them from first 
principles. 

3. Nevertheless, it may be possible to surpass the 
3% level when dealing with relative measure- 
ments of a very homogeneous data set obtained 
over a relatively short period of time and under 
similar conditions. 

4. We have derived a transformation between the 
corrected IUE data and the HST FOS absolute 
flux scale. The magnitude of the transforma- 
tions can be larger than 10% at certain wave- 
lengths. 

5. We note that much of what appears to be noise 
in NEWSIPS spectra is actually the result of 
high frequency systematics and flux calibra- 
tions. 

6. The random errors in the corrected NEWSIPS 
data are characterized in Table ^[ 

7. The pseudo trailed spectra are poorly character- 
ized by the available data. However, application 
of the large aperture corrections and flux trans- 
formation appear to reduce their systematics to 
about the 3% level. 

8. We emphasize that our results apply to NEWSIPS 
data for blue objects, and we cannot guarantee 
any broader application. 

9. Finally, a set of IDL procedures which apply the 
results of this paper to NEWSIPS low disper- 
sion spectra will be made available to the IUE 
RDAF project at the Space Telescope Science 
Institute and will be available from the authors 
on request. 
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Table 1 
Program Stars 



Star 


Spectral 
Type 


V 


B-V 


Temporal/ 
Thermal 
Standard 


Temporal/ 
Thermal 
Control 


Flux 
Standard 


Flux 
Control 


HD 60753 


B3 IV 


6.69 


-0.09 


X 






X 


HD 93521 


09.5Vn 


7.04 


-0.28 




X 




X 


BD+28°4211 


sdO 


10.52 


-0.33 


X 




X 




BD+75°325 


sdO 


9.54 


-0.37 


X 




X 




BD+33°2642 


B2 IVp 


10.84 


-0.17 




X 


X 




G191B2B 


DA 


11.78 


-0.34 




X 


X 




GD 71 


DA 


13.04 


-0.24 








X 


GD 153 


DA 


13.42 


-0.25 








X 


HZ 43 


DA 


12.86 


-0.10 








X 


HZ 44 


DA 


11.71 


-0.27 








X 



14 



Table 2 
Spectra 



Star 


Camera 


Observation 


tmin 


imax 


Number of Spectra 






Mode 


(sec) 


(sec) 


Used in Analysis 


HD 60753 


LWP 


lg apt 


4 


8 


258 






trail 


10 


50 


212 






sm apt 


o 
6 


OA 
20 


AO 

9(5 






p-trail 






7 


HD 60753 


LWR 


lg apt 


4 


20 


79 






trail 


8 


45 


77 






sm apt 


r 




22 


52 






p-trail 






2 


TTT\ 1 > / \ *— i~\ 

HD 60753 


SWP 


lg apt 


6 


12 


321 






trail 


10 


75 


270 






sm apt 


1 n 

1U 


oU 


140 






p-trail 






1 1 

_L _L 


HD 93521 


LWP 
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Table 2 — Continued 



Star 


Camera 


Observation 


tmin 


tmax 


Number of Spectra 






Mode 


(sec) 


(sec) 


Used in Analysis 


r5D+oo 2624 


aWr 


lg apt 


151) 


OCA 

350 


17y 






trail 






3 






sm apt 


250 


500 


14 
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9 
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10 


30 
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trail 


A A 

40 


1 OA 
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C A 

50 
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OA 

oU 


7A 


I I 






p-trail 






2 


nU+IO 375 
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1 a 
11) 


6b 


77 






trail 


OA 

20 


1 AA 
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20 
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75 
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Table 3 

LWP/LWR Wavelength Corrections 



A (A) 


AA (A) 


1800 


-2.50 


1950 


-2.50 


2200 


-2.00 


2250 


-1.75 


2300 


-1.50 


2385 


-1.00 


2510 


-0.40 


2730 


0.00 


3500 


0.00 



Table 4 
Wavelength Ranges 



Camera 


Range 


LWP 


1950 - 3150 A 


LWR 


1850 - 3150 A 


SWP 


1150 - 1978 A 



Table 5 

bjv1s errors and scale factors 



Camera 


Observing 


Scaling error 


a{Obs)/a(NEWSIPS) 


a(Obs)/a(NEWSIPS) 




mode 




un-normalized 


normalized 


LWP 


Large Apt 


0.020 


1.15 


1.12 




Trailed 


0.025 


1.44 


1.24 




Small Apt 


0.304 


7.93 


2.69 


LWR 


Large Apt 


0.038 


1.11 


1.09 




Trailed 


0.036 


1.18 


1.03 




Small Apt 


0.180 


3.91 


2.01 


SWP 


Large Apt 


0.018 


1.28 


1.23 




Trailed 


0.021 


1.49 


1.30 




Small Apt 


0.264 


10.3 


2.68 
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2000 
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Fig. 1 . — Ratios of white dwarf models and observations. The top panel shows the ratio of the mean IUE NEWSIPS 
data for G191B2B divided by the Finley & Koester model used to define the NEWSIPS calibration. The middle 
panel if the ratio of the HST FOS observations of G191B2B and the Koester model used in the calibration of the 
FOS. The bottom panel is the ratio of the NEWSIPS data divided by the FOS data. 
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Fig. 2. — Time dependence of standard star fluxes in the LWR camera. This figure shows the mean flux over the 
wavelength band 2400 < A < 2800 A from normalized by its mean value for 3 IUE standard stars; HD 60753, 
BD+28°4211 and BD+75°375. The mean flux was determined for each standard star observation then divided by 
the sample mean for that star. The 3 sets of relative fluxes were then overploted, with crosses for HD 60753, filled 
circles for BD+28°4211, and open circles for BD+75°375. 
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Fig. 3. — This figure demonstrates the presence of high frequency (in wavelength) time dependent systematics in 
the IUE NEWSIPS data (see 
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Fig. 4a. — Normalized SWP fluxes of IUE standard stars binned over the wavelength band 1250 < A < 1350. The 
upper plot shows the data arranged chronologically, and the lower plot shows them arranged in order of ascending 
THDA value. 
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Fig. 4b. — Same as [la], for LWP data binned over the wavelength range 2250 < A < 2300. In this case, adjoining 
points were averaged to reduce the larger noise level due to the smaller wavelength band used to demonstrate the 
effect. 
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Fig. 5a. — Examples of fits of the model to the time dependent systematics of the LWP data. The fits shown are 
for THDA = THDA , with dashed curves indicating ±0.75 of the full range in THDA. The upper left panel shows 
the fit to the relative scale factors, and the next 5 panels show fits to the data at selected wavelengths (listed on the 
figure). The wavelengths are the means of 3 IUE adjoining wavelength points to reduce the overall scatter. Each 
panel also lists the standard deviation of the 3 channel mean fit or the scale factors, a. Data from the different 
stars are keyed as follows: HD 60753 - open circles, BD+28°4211 - filled circles, BD+75°375 - filled triangles. 
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Fig. 5b. — Same as Figure pa, for LWR data. 
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Fig. 5c. — Same as Figure paL for SWP data. 



23 




2000 2200 2400 2600 2800 3000 

Wavelength 



Fig. 6a. — Ratios of fully corrected mean LWP large aperture spectra divided by FOS spectra for the 4 standard 
stars used to derive the absolute flux transformation (see Table p. The IUE wavelength scale was adjusted to agree 
with the FOS scale prior to the division. Two point binning was applied to the data for display. The individual 
ratios are depicted by different curves: solid for BD+28°4211, dotted for BD+75°375, dashed for BD+33°2642, 
and dot-dashed for G191B2B. The standard deviation of the weighted mean of the 4 curves is shown at the bottom. 
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Fig. 6b. — Same as Figure |ja] for LWR data. G191B2B is excluded due to a paucity of data. 




1200 1400 1600 1800 2000 

Wavelength 



Fig. 6c. — Same as Figure ^3a| for the SWP data. 
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Fig. 7. — Ratios of means of spectra obtained early in the mission to means of spectra taken late in the mission 
for BD+33°2642. Each set of ratios is labeled by the camera used to derive them. The dotted curves are ratios of 
uncorrected NEWSIPS data and the solid curves are ratios of NEWSIPS data corrected for systematics. The SWP 
and LWP ratios are 50 spectra means and the LWR ratios are 20 spectra means. The mean time of the spectra are 
given on the ordinate labels. 
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Fig. 8a. — Comparison of NEWSIPS and fully corrected NEWSIPS LWP data. Each pair of plots has the original, 
uncorrected NEWSIPS fluxes above the fully corrected data transformed to the FOS flux scale. 
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Fig. 9a. — Comparison of FOS data (thick curve) and fully corrected NEWSIPS LWP data (thin curve) for stars 
in common. The FOS data have been degraded to the IUE resolution. 
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Fig. 10. — Comparison of S/N values (see §0) as a function of the number of spectra averaged to make the mean. 
Large aperture spectra of HD 93521 were used in the calculations. Results determined from NEWSIPS data are 
indicated as dotted curves, fully corrected NEWSIPS results are solid curves and the theoretical, systematic-free 
limits are dashed curves. The relevant camera is indicated in each plot. 
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Fig. 11. — Ratio of the observed standard deviations, <r(Obs), to errors derived from the NEWSIPS error models, 
a(NEWSIPS), for the mean of the 3 standard stars; HD 60753, BD+28°4211 and BD+75°375. There is one panel 
for each camera-observing mode combination, and there are 2 curves for each star. The solid curves are for a(Obs) 
derived from fully corrected, unsealed observations and the dotted curve is for a(0bs) derived from fully corrected 
spectra which have been rescaled to a common mean over a fixed wavelength band. 
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Fig. 12,— Exposure level systematics in BD+28°4211 (full circles), HD 60753 (triangles), and BD+75°375 (open 
circles). The plot shows mean NEWSIPS flux values over the wavelength region 2350 < A < 2400 A for each 
star divided by the mean for all exposures with linearized flux numbers (FN) in the range 200 < FN < 400. 
Fluxes derived from saturated pixels and pixels from an extrapolation of the ITF were not included. This figure 
demonstrates that LWP fluxes derived from exposures with large FN values are systematically larger than average, 
indicating a problem with the LWP ITF (see, §0. 
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Fig. 13.— Ratios of long and short exposures for BD+28°4211 (full curve), HD 60753 (dotted), and BD+75°375 
(dashed). The plot shows the ratio of mean NEWSIPS fluxes for spectra selected to have FN values in the range 
200 < FN < 400 over the wavelength band 2650 < A < 2700 A divided by spectra selected to have FN values in 
the range 200 < FN < 400 over the same band. Data from saturated pixels and pixels using an extrapolated ITF 
were excluded. 
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